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CONVERSION FACTORS, NON-SI TO SI (METRIC) 
UNITS OF MEASUREMENT 



Non-SI units of measurement used in this report can be converted to SI (metric) 
units as follows: 



Multiply 



cubic feet per second 

feet 

feet per second 

miles (US statute) 

tons (2,000 lb, mass) 



By 



To Obtain 



0.02831685 

0.3048 
0,3028 
1.609344 
907,1847 



cubic metres per second 

metres 

metres per second 

kilometres 

kilograms 



A MATHEMATICAL STUDY OF THE IMPACT ON SALINITY INTRUSION OF DEEPENING 
THE LOWER MISSISSIPPI RIVER NAVIGATION CHANNEL 



PART It INTRODUCTION 
Backg round 

1. The US Army Engineer District, New Orleans (LMN), has been evaluating 
the feasibility of deeper draft access to the ports of New Orleans and Baton 
Rouge, Louisiana. One of the factors being considered in the evaluation is 
the impact of the deeper channel on salinity intrusion in the Lower Mississippi 
River. The importance of such an evaluation is because the city of New Orleans 
relies primarily upon the Mississippi River for its supply of fresh water. 

Approach 

2. An analysis of this factor was accomplished several years ago by the 
LMN and employed basically a steady-flow approach that assumed a channel of 
constant width and depth. The Hydraulics Laboratory, US Army Engineer Water- 
ways Experiment Station CWES), was asked if other approaches were available 
that could provide additional insight into the extent and duration of salinity 
intrusion with different low-flow hydrographs and channel sizes. Hydraulic 
Analysis Division personnel suggested a dynamic analysis of the salinity in- 
trusion phenomena using actual channel cross sections. This approach involved 
adapting an existing two-dimensional breadth-averaged hydrodynamic model to 
the Lower Mississippi River. The model (Laterally Averaged Estuarine Model 
CLAEM)) was developed by Edinger and Buchak (1981) and is the estuarine adap- 
tation of a Laterally Averaged Reservoir Model (LARM) developed for CE use 
through contract work by Edinger and Buchak (1983) and in-house work at WES. 

3. Results from on initial feasibility level analysis of the problem 
were presented by Boyd and Johnson (1982) to LMN in a letter report. In that 
report it was recommended that as the project moved into the General Design 
Memorandum phase, the analysis be repeated to develop predictions which re- 
flected the latest river geometry and a more refined model verification, 

4. To provide new river geometry data at key river crossings, LMN 



personnel conducted hydrographic surveys in October 1982. During the analysis 
of these data, it was discovered that several major river crossings below New 
Orleans had not been accurately represented in the preliminary analysis. 
Since such crossings can play a major role in retarding the intrusion of the 
salt wedge, this reemphasized the need to reanalyze the impact of a deeper 
channel on salinity intrusion using recent geometry data at the crossings. 
Also, in the initial analysis, only a very limited study of the sensitivity of 
model results to diffusion coefficients and downstream boundary conditions on 
salinity and the water-surface elevation was attempted. In the study reported 
herein, a more detailed analysis of these factors is presented. 

Purpose 

5. A major reason for conducting the numerical modeling study was not 
only to evaluate the impact of a deeper channel on salinity intrusion using 
historical riverflows but also to use the numerical model to determine the 
impact on the wedge intrusion of increasing the height of natural river cross- 
ings, i.e., creating a sill in the river. To address the question of the 
stability of such a sill composed of natural sediment, the HEC-6 computer pro- 
gram "Scour and Deposition in Rivers and Reservoirs," HEC (1977), was applied. 



PART III DESCRIPTION OF NUMERICAL FLOW - SALINITY MODEL 

6. Under the Environmental and Water Quality Operation Study (EWQOS) 
program of the US Army Corps of Engineers, developmental work has been con- 
ducted for several years on a two-dimensional (2D) laterally averaged, free 
surface, variable density, and heat conducting model for use in simulating 
flows in thermally stratified reservoirs. This effort, which extended the 
earlier work funded by the US Army Engineer Division, Ohio River (Edinger and 
Buchak 1979) , has resulted in a numerically efficient model that is known as 
LARM. Under a contract with the US Army Engineer District, Savannah, LARM has 
been modified for use in computing stratified flows in estuaries as a result 
of both salinity and thermal effects. This model is known as LAEM, Edinger 
and Buchak (1981), and was the model applied in the present study. 

Governing Equations 

7. The basic set of equations that are solved in LAEM are statements 
of the conservation of mass and momentum of the flow field plus the conserva- 
tion of the heat and salt content of the water body. The governing equations 
are developed by first performing a temporal averaging of the three-dimensional 
equations for laminar flow. Boussinesq's eddy coefficient concept is then em- 
ployed to account for the effect of turbulence in the flow field. Next, the 
time-averaged equations are averaged over the estuary width and finally over 

an individual vertical layer with boundaries at z = k + 1/2 and z = k - 1/2 
to yield the following equations that are solved in LAEM: 

longitudinal (x-direction momentum) 



£ CUBh) + £ (u 2 Bh) + (u b w b b) k+1/2 - Cu b w b b) k _ 1/2 
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+ £ £ (PBh) - A x 72 < UBh) + < T z b >k+l/2 - ( ^ b) k-l/2 = ° < U 



with 



T = C*p„/pW cos 6 (surface) 
z d a 



= A 3U/3z (interlayer) 
z 

- gUJU|/c 2 (bottom) 



internal continuity 

Cw b b) k-l/2 = (w b b) k+l/2 + h (UBh) " qBh C2) 

total depth continuity 

W ' I h < UBb > = J ^ ( 3) 

k k 

vertical ( z-direction? momentum 
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3^=Pg (4) 



heat balance 
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salinity balance 
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equation of state 

p = TlOOO P /(LA + 0.698 P fl (7) 

*■ o oJ 



where 

P = 5890 + 38T - Q.375T 2 + 3S 
o 

LA = 1779.5 + 11.25T - 0.0745T 2 (8) 

- (3.8 + 0.01T)S 

Variables in the equations above are defined below: 

2 

A x-direction momentum dispersion coefficient, m /s 
x 

2 
A z-direction momentum dispersion coefficient, m /s 

b Estuary or river width, m 

B Laterally averaged width integrated over h , m 

1/2 

c Chezy resistance coefficient, m /s 

C* Resistance coefficient associated with wind 

2 
D x-direction temperature and salinity dispersion coefficient, m /s 

2 
D z-direction temperature and salinity dispersion coefficient, m /s 

2 
g Acceleration due to gravity, m/s 

h Horizontal layer thickness, m 

H Source strength for heat balance, Cm s 
n 

k Integer layer number, positive downward 

}' Pressure, N/m 

q Tributary inflow or withdrawal, m /s 

S Laterally averaged salinity integrated over h , ppt 

t Time (s) 

T Laterally averaged temperature integrated over h , C 

u, x-direction, laterally averaged velocity, m/s 
b 

II x-direction, laterally averaged velocity integrated over h , m/s 

8 



3 

V Cell volume (B • h • Ax) , m 

W Wind speed, m/s 
a 

w, z -direction, laterally averaged velocity, m/s 
b 

x and z Cartesian coordinates: x is along the estuary center line at 
the water surface, positive to the right and z is positive 
downward from the x-axis, m 

Ax Longitudinal spatial step, m 

At Timestep, sec 

5 Surface elevation, m 

3 

p Density, kg/ra 

3 
p Air density, kg/m 

a- 

t Tangential stress in positive x-direction divided by the water 

Z 2 2 

density, m**/s 

$ Wind direction, rad 

8. Basic assumptions in addition to the reduced dimensionality are that 
the Boussinesq approximation (p is constant except where multiplied by 
the acceleration of gravity) is applicable and that vertical accelerations 
are negligibie so that the pressure can be considered hydrostatic. In addi- 
tion, the concept of eddy coefficients is utilized to represent the effect of 
both time averaging, as previously noted, and spatial averaging of the equa- 
tions. The horizontal dispersion coefficients, A and D , are assumed to 

xx 

be constant, whereas the vertical dispersion coefficients, A and D z , are 
dependent upon the stratification as reflected by the local Richardson number, 

A = A (1 + 3.33R)" 3/2 
z z 
o 

D - D (1 + 1GR.)" 1 ' 2 (9) 

z z 1 

o 



(R i ) i.e 



where 



(0 



R. - "f 



and A and D are the vertical coefficients for no stratification. 

z o z o 
Due to the hydrostatic pressure assumption, unstable stratification cannot 



be modeled in a convective fashion and thus is handled in a diffusive manner 

2 
by increasing D to its stability limit of h /2At , where it is the 

computation time-step. 

9. The laterally averaged horizontal pressure gradient in the longitu- 
dinal momentum equation contains the density driving force. Using the expres- 
sion for the hydrostatic pressure, the horizontal pressure gradient can be 
divided into its two components of the barotropic (surface slope) gradient 
and the baroclinic (density) gradient to yield: 



£♦■/&) 



If " "8" & 



Numerical Solution Scheme 

10. Finite difference techniques are employed to solve the governing ' 
Equations 1-7. The particular scheme employed is structured such that the 
water-surface elevations are computed implicitly. Using the new water-surface 
elevations, the x-component of the flow velocity is then explicitly computed 
from the longitudinal momentum equation. As in other hydrostatic models, the 
vertical component of the velocity is computed from the continuity equation 
which is reduced to the incompressibility condition as a result of the 
Boussinesq approximation. The solution begins at the bottom and progresses 

up the column of layers. With the flow field computed, the water temperature 
and salt concentration are then computed from their respective transport 
equations in an implicit fashion. A detailed discussion of the solution 
scheme can be found in the report by Edinger and Buchak (1979) . 

11. The major advantage of the solution scheme employed in LAEM is 

that the extremely restrictive stability criterion based upon the speed 

of the free surface gravity wave, i.e. At < (Ax/VgH ) , where H is 

max max 

the maximum water depth, is removed. However, since the convective terms 
and density gradient in the longitudinal momentum equation, as well as the 
vertical diffusion term in the temperature and salt transport equations, are 
lagged in an explicit fashion the following stability criteria still remain: 

10 



i«<* 



At < 



Ax 



^T 



(12) 



and 



At < 



2D 



where h is the layer thickness, Ap is the density difference of the fresh 

and saline water, and h is the height of the density flow. 

w 

12. The finite difference grid employed in LAEM is the marker and cell 
grid employed by Leendertse (1967) and others in 2-D vertically averaged 
models. As shown below, the velocity components are defined on the faces of 
a cell, with the water temperature, salinity, and corresponding density defined 
at the center 

COMPUTATIONAL GRID 



k - 1/2 



k + 1/2 




Values of variables needed at locations where they are not defined are ob- 



tained by averaging, e.g. 



w =1/2 (w, + w, ) . 

cell center k-1/2 k+1/2 

Data Requirements 



13. The major data input required by LAEM is the geometry data des- 
cribing the system. At the center of each computational cell the width of 



11 



the estuary or river must be prescribed. For the study described herein, 
these were obtained from hydrographic survey data furnished by the Lower 
Mississippi Valley Division (LMVD) and LMN. Additional data required are the 
boundary conditions that drive the internal flow field. At the river bound- 
ary, a discharge hydrograph must be prescribed, whereas at the ocean bound- 
ary the tide must be input. In addition, at inflow boundaries the temperature 
and salinity must be prescribed. In the present study, the water temperature 
was assumed constant and the surface heat exchange was set to zero. However, 
for problems in which thermal effects are considered, shortwave solar radia- 
tion, air temperature, dew-point temperature, and wind speed must be known 
in order to compute the coefficients required in the computation of the rate 
of surface heat exchange. 



12 



PART III: SYSTEM DESCRIPTION AND SCHEMATIZATION 
System Schematization 

14. Figure 1 is a map showing the Mississippi River from Baton Rouge 
at mile 228,4 AHP* to the Gulf of Mexico. This segment of the Mississippi 
River, which would be deepened and/or widened where necessary to provide a 
750-ft-wide** by 45-, 50- or 55-ft-deep navigation channel (with overdepth 
dredging), was schematized in the modeling effort (Figure 2). At the Gulf where 
flow exits through three major distributaries (Southwest Pass, South Pass, and 
Pass a Loutre) and several smaller ones, it was assumed that the system could 

be satisfactorily reproduced by treating Southwest Pass as the dominant dis- 
tributary and representing outflow through the other distributaries as lateral 
outflows at the appropriate locations. This assumption also implies that 
wedge intrusion into the system is dominated by saline water moving up South- 
west Pass. These assumptions were discussed with WES and LMN personnel, who 
have extensive experience with the passes area, and appear reasonable for the 
existing situation. They will be even more appropriate under future condi- 
tions as South Pass continues to shoal from the earlier maintained navigation 
depth of 30 ft to a controlling depth of about 17 ft, similar to the control- 
ling depth that now exists in Pass a Loutre. The dredging of a deeper naviga- 
tion channel would also tend to make the assumptions more accurate. 

Flow Distribution in Distributaries 

15. Using the approach described above, it was necessary to prescribe 
the distribution of flow through the different distributaries. Discussions 
with WES and LMN personnel provided general guidance concerning the appropri- 
ate distribution. A one-dimensional unsteady flow model, System 11, devel- 
oped by the Danish Hydraulic Institute, has been applied to the passes area.t 
This work, which used the most accurate geometric data available to describe 



* AHP - Above Head of Passes. 
** A table of factors for converting non-SI units of measurement to SI (metric) 

units is presented on page 3, 
t Personal communication, R. Athow (1981), Hydraulics Laboratory, US Army 

Engineer Waterways Experiment Station, Vicksburg, Miss. 
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Figure 1. Location map of study reach 
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Figure 2. Schematic of Lower Mississippi River 
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flow distribution 



the complex distributary system, provided estimates of outflow in each 
distributary that were used in the 1982 study by Boyd and Johnson. However, 
it was decided that the flow distribution more recently provided by LMN would 
be used in the current work by making the assumption that the distribution 
provided would be reasonably appropriate for the discharge range of interest 
(roughly 100,000 to 300,000 cfs). The distribution specified was: 

Distributary Approximate River Mile Percent of Flow 
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West Bank 


3 


2.3 


South Pass and 






Pass a Loutre 





45.4 



As illustrated in Figure 3, this leaves 30.2 percent of the riverflow to be 
discharged through Southwest Pass as opposed to 36.8 percent in the earlier 
study. 

Channel Geometry 

16. Hydrographic survey information obtained during the period 1973-1975 
was used to develop channel sections for the feasibility study model. Cross- 
section data at intervals from about 1/2 mile to 3 miles were furnished by LMVD 
and/or LMN. To conserve computer time, the vertical grid used in the model 
ended at a depth of 100 ft. Plots of the individual cross sections were in- 
spected, and where appropriate, adjustments were made in channel widths above 
-100 ft to roughly conserve volume. These adjusted cross sections were used in 
the GEDA computer code ( HEC 1976) to develop the needed channel elevation-width 
tables at 1-mile intervals along the river between the Gulf and Baton Rouge. 

17. Since major crossings in the river can have a significant impact 
on the salinity intrusion, it was considered necessary in the reanalysis to 
have an accurate representation of these crossings. A major question was 
whether or not the crossings reflected in the 1973-75 survey had experienced 
significant changes in the past 8 to 10 years. To answer this question, LMN 
conducted a detailed hydrographic survey of the five major crossings that 
would be likely to influence the intrusion of a salt wedge during low-flow 
periods. Cross-sectional plots of these crossings are presented in Figures 4-8, 
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Figure 4. Crossing at mile 115.9 
from 1982 and 1973-75 data 



a 
> 
a 

2 



Z 

o 
p 

ID 



D 


- 

1 








-10 


- 








-20 










-30 




>~-1 949-52 


1 


-40 






If 
1 
if 


-so 

-60 


_ 


T 1 


/^1982 

/ /J 
/ / 
* s 

-* y~— ^ 

1 I 1 1 



500 



1000 1500 2000 
WIDTH. FT 



2500 3000 



Figure 5. Crossing at mile 89.7 
from 1982 data and 1949-52 data 
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Figure 6. Crossing at mile 79.2 
from 1982 data 
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Figure 7. Crossing at mile 63.4 
from 1982 data 
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Figure 8. Crossing at mile 47.3 from 
1982, 1973-75, and 1949-52 data 



with Figures 4, 5, and 8 illustrating 
the relatively strong stability of the 
crossings through comparisons of 1982 
sections with 1973-75 and/or 1949-52 
sections. 

18, Figure 9 presents a schematic 
of the adjusted channel bottom profile 
which reflects the lowest active cell 
in the computation grid at each point 
along the channel. Of course, channel 
widths are highly variable with eleva- 
tion and along the channel. As men- 
tioned earlier (paragraph 4) , it was 
discovered during the analysis of the 
new survey data that the feasibility 
level study representation of some 
crossings did not accurately reflect 
the 1973-75 hydrographic data. Fig- 



ure 9 shows the additional crossings that have been added to the model since 
the earlier study. With a longitudinal spatial step, Ax , of 2 miles and a 
vertical spacing of 5 ft, the computational grid is composed of (27 x 126) 
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c ells. In the earlier study, initial runs were made with Ax = 1 mile. How- 
ever, computer experimentation revealed essentially no difference in either the 
flow or salinity computations when Ax was increased to 2 miles; therefore 
the larger value has been used in subsequent computations to save on computer 
costs. 

19. Representation of the proposed deeper channels required deepening in 
Southwest Pass (indicated by the dashed line in Figure 9 for a 55-ft channel) 
and widening of a number of sections along the channel to provide the desired 
750-width. 
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PART IV: MODEL VERIFICATION 



Available Field Data 



20. Since the occurrence of the 1973 flood on the Lower Mississippi 

River there has been no dredging of South Pass, with the channel subsequently 
shoaling to its current depth of about 17 ft. Recall that this resulted in 
the basic assumption that the salt wedge is currently primarily driven by 
saline waters moving through Southwest Pass which resulted In the assumption 
that the study reach could be modeled as a single channel extending from Baton 
Rouge through Southwest Pass to the Gulf. With this assumption, field data 
on wedge movement collected before 1973 were of limited use in verifying 
the model. Therefore the verification of the model was forced to rely upon 
a limited set of data collected during a relatively low-flow period in 1981. 
This set of data consists of the time-history of the tip of the wedge as well 
as vertical profiles of the salinity behind the tip of the wedge in January 
1981. In addition, it has been generally observed that the tip of the wedge 
(assumed to be 9 ppt salinity or 5 ppt chloride) intrudes to the vicinity of 
the Head of Passes for a riverflow of about 300,000 cfs. Initially, an effort 
was made to ensure that the computed water-surface profiles were in reasonable 
agreement with field data. Stage information at Baton Rouge, New Orleans, and 
Venice were available for the 1980-1981 hydrograph, whereas estimates of 
appropriate stages for a steady riverflow of 300,000 cfs were developed using 
stages observed during September 1975 when the riverflow fluctuated around 
300,000 cfs. 

S ensit i v it y Runs 

Influence of additional crossings 

21. Figure 9 shows that the major crossing at about mile 47 was not 
well represented in the earlier study. Therefore the verification of the model 
from the earlier study, which centered around forcing the tip of the wedge to 
stabilize at about mile 60 for a steady riverflow of 150,000 cfs, was unlikely 
to be appropriate now. Figure 10 illustrates that this was indeed true since 
the crossing at mile 47 completely stopped the wedge for the 1980-1981 river- 

f lows . 
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Figure 10. Effect of crossing at mile 47.3 
and a variable salinity condition 

Influence of salinity profiles used 
a s boundary conditions at the Gulf 

22. In the earlier study, the saLinity profile used as the boundary con- 
dition at the Gulf was constructed from previous physical model tests. The 
Profile corresponding to a riverflow of about 125,000 cfs was selected to be 
used over the complete range of flows encountered in the model runs (about 
90,000 to 300,000 cfs). In the present study, salinity profiles corresponding 
to river flows of 125,000 and 350,000 cfs were employed (Figure 11) with linear 
interpolation providing the boundary salinity at intermediate flows. For 
r iverflows below 125,000 cfs, the profile corresponding to 125,000 cfs was 
Prescribed since no additional data were known. As shown in Figure 10, the 
time required for the wedge to reach the crossing at mile 47 is significantly 

-7,5 




20 15 

SALINITY. PPT 

Figure 11. Variable salinity boundary 
condition at the Gulf 



21 



influenced when employing the variable boundary salinity condition for the 
1980-1981 flows. 

23. It should be noted that since the advective terms in the salinity 
equation are represented by one-sided differences based upon the direction of 
flow, only those values of the prescribed salinity profile lying in the lower 
cells where the denser saline water moves up the river enter into the computa- 
tions. Therefore, in both a mathematical and physical sense, the downstream 
salinity boundary condition enters the computations properly. 

Influence of vertical eddy coefficients 

24. Both the x-momentum equation and the salinity transport equation 
contain eddy coefficients in both the x~ and z-directions. Two orders of 
magnitude change In x-direction coefficients result in very little change In 
the computed flow and salinity fields. However, the model Is more sensitive 
to variations in the vertical coefficients. As shown In Figure 12, by de- 
creasing the vertical coefficients <D - salinity equation, A - momentum 
equation) by a factor of 6, the wedge can be farced to move over the crossing 
at mile 47 an additional 8 miles for the 1980-1981 flows. Figure 13 illus- 
trates the Influence of the vertical eddy coefficients on the salinity In- 
terface corresponding to 9 ppt. These results show that the more dominant 
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of the two coefficients in influencing salinity intrusion is the vertical dif- 
fusion coefficient, D , in the salinity transport equation. Relatively 
large variations in A seem to have little influence on the computed flow 
field and thus little influence on the salinity field. Edinger and Buchak 
(1981) arrived at a similar conclusion in sensitivity studies using data from 
the Potomac River. 
Influence of Ghezy coefficient 

25. The major coefficient to be varied to match a computed flow field 
with observed data is the Chezy coefficient, C . Figure 14 illustrates 

that the value of C has a major influence on the water-surface profile com- 

Z 1/2 1/2 

Puted . Varying C from 45 m /sec to 100 m /sec results in a decrease of 

about 6 ft in the water surface at Baton Rouge and about 3 ft at New Orleans 
for a steady flow of 300,000 cfs and a Gulf elevation of 0.0 ft. Since the 
flow and salinity fields are coupled through the longitudinal pressure gra- 
dient, major changes in the flow field result in significant changes in the 
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Figure 14. Effect of C on water-surface profile 
for 300,000-cfs steady flow 

salinity computation. Figure 15 shows that for a steady flow of 300,000 cfs, 
the salt wedge only moves about 2 miles up Southwest Pass if C =45 through- 
out the system, whereas the wedge stabilizes at a distance of more than 

26 miles from the Gulf for C = 100. 

z 

Influence of Gulf elevation 

26. The 1980-1981 hydrograph (Figure 16) was used in runs of the final 
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wedge interface (9 ppt) for 300,000-cfs 

steady flow and a variable salinity Figure 16. 1980-81 discharge used 

boundary condition as boundary condition at Baton Rouge 

verified model to investigate the influence of the Gulf elevation boundary 
condition on the salinity intrusion. In the earlier application of LAEM, the 
Gulf elevation was set to be constant at 0.0 ft NGVD after runs with a sinu- 
soidal tide revealed no major impact on movement of the salt wedge. Therefore 
the initial verification of the model with the additional crossing data em- 
ployed a 0.0 ft Gulf elevation also. However, after talks with WES personnel 
revealed that data they had analyzed seem to indicate that long term "setups" 
or "setdowns" occur in the Gulf, it was decided to make additional runs to 
determine the influence on the wedge movement of using recorded elevations in 
the Gulf rather than assuming a zero elevation. To accomplish this, the re- 
corded elevations presented in Figure 17 at the Grand Isle gage were employed 
as the downstream boundary condition, with the discharges shown in Figure 16 
as the upstream boundary condition. The location of the Grand Isle gage is 
some 30 to 40 miles west of the end of Southwest Pass, as illustrated in 
Figure 1. 

27. Comparisons of computed elevations, using both a zero Gulf eleva- 
tion and the Grand Isle elevations as downstream boundary conditions, versus 
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recorded data at Baton Rouge, New Orleans, 
and Venice are presented in Figures 18-20. 
In addition, a comparison of the water- 
surface profiles computed for the two 
cases at the end of the hydrograph when 
the riverflow has returned to about 
300,000 cfs is presented in Figure 21. 
Even though the use of a zero Gulf eleva- 
tion does not reflect the oscillations 
that are obviously occurring, in a mean 
sense the comparisons of recorded and com- 
puted elevations at Baton Rouge and New 
Orleans are in reasonable agreement while 
the computed elevations at Venice are 
about 0.5 ft too low. However, this is 
not the case at New Orleans and Venice 
when the recorded Grand Isle elevations 
Were employed as the boundary condition 
at the Gulf. Making adjustments in the 

Chezy coefficient to force the water-surface elevations at New Orleans and 
Venice to match the recorded data would result in little hope of ever forcing 
the movement of the salt wedge to match the field data. Therefore one is faced 
with the question of whether or not to "believe" the Grand Isle data. The os- 
cillations in the water surface, which may occur over several days, obviously 
are real. Therefore the basic problem seems to revolve around determining the 
proper zero of the gage, i.e., it would seem that the curve in Figure 17 should 
be shifted upward to yield a mean elevation of approximately zero. If it is 
accepted that the mean Gulf elevation should be 0.0 ft NGVD, the next question 
to be addressed is the influence of the relatively long-term oscillations ob- 
served in the Grand Isle data on the wedge movement. This question was an- 
swered by imposing a sinusoidal oscillation with an amplitude of 2.0 ft and a 
Period of 10 days at the Gulf. The effect of such long-period oscillations 
on the wedge movement is illustrated in Figure 22. As can be seen, the long- 
term effect over the hydrograph is negligible and had little influence on the 
extent and duration of the salt intrusion. Therefore it seems reasonable to 
conclude that the mean water-surface level of the Gulf rather than the tidal or 



Figure 17. Recorded water- 
surface elevations at Grand 
Isle 
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Figure 18. Comparison of computed 
and recorded water-surface eleva- 
tions at Baton Rouge 
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and recorded water-surface eleva- 
tions at Venice 
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longer-period oscillations exerts the major influence on the movement of the 
s alt wedge. Based upon these results, it seems reasonable to conclude that 
the use of a zero Gulf elevation as the downstream boundary condition is the 
a ppropriate condition to use in both the verification run and the series of 
historical hydrographs to be discussed later. 
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Verification Results 

28. As illustrated in Figure 10, after placing the additional crossings 
In the system, the model as calibrated from the earlier study could not repro- 
duce the 1980-1981 field data on the wedge movement. In the earlier study, 
the final Chezy adjustment used a value of 100 in Southwest Pass and a value 
of 45 above Head of Passes. With the new crossing at mile 47, a problem in 
the verification exercise developed. Field data indicate that the wedge moves 
past the crossing and essentially stabilizes at about mile 61. Therefore ad- 
justments in the Chezy coefficient and/or the vertical diffusion coefficient, 

D , are required to force the wedge to overtop the crossing at mile 47 

o 

fairly quickly; but the adjustments must be such that once the crossing is 

overtopped the wedge stops within the next 10 to 15 miles. 

29. Many different combinations of the roughness and diffusion coeffi- 
cients were tried. The final values of A and D selected were 

-4 2. -4 2, Z o Z ° 

1.5 x 10 m /sec and 0.5 x 10 m /sec, respectively. The model default value 

-4 2 
for both is 2.0 x 10 m /sec, which is representative of values in the Po- 
tomac River. 

30. With the spatial variation of C shown in Figure 9, the values of 
A and D noted above, the variable salinity boundary condition given in 
Figure 11, and a zero Gulf elevation as the downstream elevation boundary con- 
dition, the final verification results on wedge movement for the 1980-1981 
hydrograph are presented in Figure 23. A physical explanation for the large 
difference between Chezy values above and below Head of Passes is not apparent. 
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It could indicate that the assumed flow diversions in the Passes area are not 
entirely correct. It might also be noted that a simplified form of the x- 
raomentum equation is solved at the downstream boundary, which, as noted by 
Edinger and Buchak (1981), may result in the need for a nonphysical value of 

the Chezy coefficient near the downstream boundary. Using the expression 

1/6 
n = h /c to convert the Chezy coefficient to the more commonly employed 

Manning roughness coefficient in river studies yields corresponding n-values 

°f 0.015 in Southwest Pass and about 0.035 in the river. 

31. Figures 18-21 (as discussed earlier) show the comparison of computed 
and recorded elevations at Baton Rouge , New Orleans, and Venice plus the water- 
surface profile for a riverflow of approximately 300,000 cfs at the end of the 
1980-1981 hydrograph. Remember that the Gulf elevation is 0.0 ft. 

32. Computed and observed vertical salinity profiles are plotted in 
F *gure 24. The field data plotted are from measurements on 28 January 1981 
at a location 5 miles behind the tip of the wedge, which is located at about 
mile 59. The computed profile presented is on the same day, i.e., 28 January 
-1981 , but at a distance of 8 miles behind the wedge tip. As can be seen, 
these results compare reasonably well with the field data. As a result of 
numerical diffusion in the model (exists in some degree in virtually all 
numerical models), the computed interface is not as sharp as the field data 
imply. However, the computed and recorded concentration of approximately 

l *G ppt that occurs at a depth of about 50 ft below the water surface compares 
We ll. Additional discussion on the vertical structure of the wedge is provided 
in PART V. 
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33. The model as adjusted responded properly to boundary condition 
changes, and in general, the agreement between computed results and 1980-1981 
field observations used in the verification was quite good. A set of field 
data for a much greater intrusion of the wedge would have been preferable; 
however, as previously noted, these were not available for the system as it 
exists today. 
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PART V: ARRESTED WEDGE ANALYSIS 

34. Results from computations made with a constant riverflow of 
150,000 cfs are used to analyze the characteristics of the arrested saline 
wedge formed. These are then compared with experimental results obtained by 
Keulegan ( 1949) in laboratory studies. The vertical distribution of salinity 
and the longitudinal component of the velocity, as well as the position of the 
arrested wedge interface and the inter facial shear stress, are considered. 

Distribution of Salinities 

35. For a steady flow of 150,000 cfs, the toe of the arrested saline 
wedge stabilizes at mile 61 above the Head of Passes. The computed vertical 
distributions of salinity at miles 8, 34, and 54 are presented in Figure 25. 
The common procedure is to plot the salinity as a ratio of the salinity at a 
particular depth and the bottom salinity, S . The interface surface is then 
taken where S/S equals 0.5. The freshwater depth, h. , is the distance 
from the surface to that point. Figure 25 shows that the freshwater depths at 
miles 8, 34, and 54 are 38.0, 61,3, and 68.2, respectively. 

36. An effective interface layer thickness, z , is introduced through 
the relation below 



(#) 



2 C \ s sr L ' > <13) 

max 



where z is positive downward. For the three locations plotted in Figure 25, 
the values of z are determined to be 18.52, 20.0, and 13.5, respectively. 
The ratio, h /z , is then computed to range from 2.05 to 5.05, with the aver- 
a ge value being 3.39. Keulegan determined a value of 14.4 in his experiments 
with fresh water flowing over a stagnant saltwater layer. From Figure 26, 
which is a plot of field data collected in 1981 at mile 22 and mile 55, the 
values of h./z are determined to be 4.02 and 6.90, respectively. Therefore, 
although Keulegan' s experiments imply that the numerical results at the 
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interface are much too diffuse, the limited field data available agree relative- 
ly closely with the computed results. 

37. The use of z leads to a dimension less representation of the 

P 

salinity. In Figure 27, S/S is plotted against the ratio, z'/z , where 
z' is the vertical distance from the interface line, for miles 8, 34, and 54. 
As can be seen, similarity is preserved rather well and the line drawn through 
the computed data points is the function 



(' " ta " h v) 



S/S q - 1/2 [l - tanh ^-\ (14) 



which is the same function determined by Keulegan (1949). 

Distribution of Velocity 

38. The vertical distribution of the horizontal component of the flow 
velocity in the three cross sections at miles 8, 34, and 54 are given in 
Figure 28, where the ratio of the velocity, U , to the surface velocity, U , 
is plotted against 2'/n. . Although there is some deviation, similarity seems 
to be somewhat preserved except for the velocity in the saline layer at mile 8. 
However, it should be noted that the saline layer at mile 8 lies in the hole 
created by vertical walls at miles 5 and 11 (see Figure 9) and thus might be 
expected to behave differently. 

39. An expression for the representative curve drawn in Figure 28 may 
be obtained as follows. Since the velocity close to the surface changes slow- 
ly, it is assumed that 

d(U/U ) 

~~-= CI - I}* <a o + a in ) (15) 

where 

n = z'/h 1 

a = value of the rate of change of velocity at n = 

a. = constant to be determined 
From Figure 28, a is determined to have a value of 2.5. Integrating, and 
using the conditions that U/U = 0.25 at n = and U/U =1.0 at r\ = 1 , 
yields 



33 




-1.0 



Figure 27. Distribution of normalized salinity about the 

wedge interface 



1.0 - 
0.8 - 

0.6 - 

0.4 

0.2 

c 


-0.2 - 
-0.4 - 

-0.6 

-08 









i 


— 






o7 


- 






o 7 

/• 


- 




















*M 






o 






~ 


/i 






_ 


r 




LEGEND 




V 




• - MILE 8 




■ \ 




* - MILE 34 




• 




o - MILE 54 




i • i i 


1 


1_ i i \ 



-0.4 -0.2 0.4 0.6 0.8 1.0 
U/U s 

Figure 28. Distribution of normalized longitudinal 
velocity component about the wedge interface 



34 



U/U = 1/4 + 5/2n - 3n 2 + 3/2n 3 - l/4n 4 (16) 

s 

Now, using the expression for U/U , the average freshwater velocity can be 
obtained from 

1 

u = /u/u s dn (17) 

o 

to yield 

U = 0.83 U (18) 

A numerical integration of the computed freshwater velocities gives 

U = 0.75 U (19) 

s 

Therefore the polynomial form for the freshwater velocity seems appropriate. 

40. Keulegan's experiments gave a value of 0.54 U for the inter- 
facial velocity, whereas Figure 28 yields a value of 0.25 U for the numer- 
ical results. The difference can probably be attributed to the much sharper 
interface in Keulegan's experiments. 

Interfacial Stress 

41, Neglecting the momentum of the saline water and sidewall friction 
effects, Keulegan gives the equation below for the interfacial shear stress, 



(20) 



where 

U = average freshwater velocity 

h. - depth of the freshwater layer 

Ap = maximum density difference between the saline and fresh water 

h_ ■ depth of the saline layer. 
Figure 29 is a plot of Ap/p versus salinity, as determined from Equation 7, 
with the water temperature set to 20 C. The first step in the determination 
of t, is the evaluation of dh„/dx . Since 
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1 + h 1 /h 2 



4* 
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Figure 29. Ratio of excess water density to fresh- 
water density as a function of salinity at a 
temperature of 20 C 



100 + 6 = h + h 



(21) 



dh ? /dx can be evaluated from 



dh 2 /dx = -dhj/dx + d6/dx 



(22) 



where 6 



dh„/dx , where 



is the water-surface elevation above mean sea level. Values of 

x has been taken as positive in the downstream direction, as 

are given at various 



well as other variables needed in the evaluation of t 
river miles in Table 1. Performing the computations in Equation 20 and aver- 
aging over the arrested wedge yields an average value of 



'avg 



= 0.00107 



(23) 



Following Keulegan (1949), the interfacial friction coefficient is defined to 



be 



f. *= S(t./pU 2 ) 
i \ i /avg 



(24) 



which yields an average value of f 



0.0086 



However, recall that the inter- 



face was defined as that point where the salinity is 1/2 of the bottom salinity 
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If the value for Ap/p used in the computation for t. had been based upon 
the interface salinity rather than the bottom salinity, the computed value of 
f. would have been decreased by a factor of about 2, i.e. the average computed 
f . would have been approximately 0.0043. 

Critical Flow at the Gulf 

42. One important computation for the arrested wedge remains to be made. 
It is generally assumed that the river mouth is a critical flow section with 
the densimetric Froude number being equal to unity for the case of highly strat- 
ified salt wedges in estuaries. In other words, the average freshwater veloc- 
ity at the mouth is 



°->pr^ 



(25) 



For the case of a constant 150,000-cfs flow, the computed freshwater depth was 
determined to be 10.39 ft and Ap/p = 0.024 . Therefore the average freshwater 
velocity at the river mouth should be 2.83 fps. The numerical results yield 
an average velocity of 2.53 fps which results in a densimetric Froude number of 
0.89. Therefore the numerical results imply that the river mouth is close to 
being a critical flow section with respect to the densimetric Froude number. 

43. The discussion above has been for the purpose of comparing LAEM 
computed results for an arrested saline wedge with experimental results from 
Keulegan and known theoretical results. In general, the results compare quite 
well and serve to increase the credibility of the unsteady results obtained 
from LAEM that are discussed in PART VI of the report. 
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PART VI: APPLICATION OF MODEL TO HISTORICAL FLOWS 
Steady-Flow Tests 

44, An analysis of the impact of a 55-ft channel on salinity intrusion 
in the Lower Mississippi River was accomplished several years ago by LMN and 
employed basically a steady-flow approach which assumed a channel of constant 
width and depth. In the earlier study by Boyd and Johnson (1982), steady-flow 
tests were run for riverflows of 300,000, 150,000, 125,000, and 91,600 cfs. 
The last three flows correspond to the low flows experienced during 1980-81, 
1963, and 1936, respectively. 

45. The 300,000-cfs tests were initiated from a no-flow condition with 
the flow being Increased over a few hours to 300,000 cfs and then held at 
300,000 cfs until the wedge position stabilized. Other steady-flow tests were 
initiated with a flow of 300,000 cfs or less, changed over a few hours to the 
desired flow, and held there until the wedge position stabilized. The posi- 
tions at which the wedge stabilized for the different channels and riverflows 
are plotted In Figure 30 and tabulated below. 

40-ft Channel 55-ft Channel 



Wedge Location 


Computation 


River Mile 


Length, days 


11 


10 


92 


30 


116 


30 


127 


50 



Tarbert Landing Wedge Location Computation 
Discharge, cfs River Mile Length, days 

300,000 5 

150,000 58 24 

125,000 93 28 

91,600 116 59 

The computed wedge intrusion is generally supportive of the curves developed 
by LMN (Figure 30), particularly at the lower flows. For the existing 40-ft 
channel and a riverflow of 91,600 cfs, the wedge was stopped by the Kenner 
Crossing at mile 116. The computation was run a total of 59 days to make sure 
that the computation would not ultimately allow the wedge to spill over and 
move farther upstream. For a proposed 55-ft channel, the wedge was stopped by 
Kenner Crossing for a discharge of 125,000 cfs. When the discharge was de- 
creased to 91,600 cfs, the wedge did move past the Kenner Crossing but its ul- 
timate position (mile 127) was greatly Influenced by the crossing. As noted, 
these results were obtained from the earlier study In which the additional 
crossings were not considered. Based upon results from the hydrograph runs 



38 



300 r 



250 



- 200 

ui 

» 

E 

< 

i 

Q ISO 
U 

5 

Z 



100 



BO 



LMN E5-FT CHANNEL 



-LMN 40-FT CHANNEL 



« COMPUTED 40-FT CHANNEL 
• COMPUTED S6-FT CHANNEL 



140 



120 



100 SO 60 40 20 

POSITION OF tVEDGE, MILES ABOVE HEAD OF PASSES 



-20 



Figure 30. Wedge position for steady-flow tests 
discussed next, in which the crossings were considered, essentially the same 
results would be obtained, 

Hydrograph Tests 

46. Hydrograph tests were run to provide insight concerning the extent 
and duration of wedge intrusion for changing riverflow conditions. LMN ex- 
pressed interest in studying eight low-flow hydrographs which represented a 
range of frequency of occurrences and shapes. These hydrographs were 1936, 
1938-39, 1947-48, 1952-53, 1953-54, 1955-56, 1968, and 1980-81. All hydro- 
graphs were initiated from a riverflow of 300,000 cfs with initial conditions 
that had been saved from a 300,000-cfs steady-state run of the most recently 
calibrated model containing the additional crossings. Transient solutions 
were then computed using a time-step of 15 tnin with recorded Tarbert Landing 
discharges, lagged by 2 days, prescribed as the upstream boundary condition at 
Baton Rouge and a zero Gulf elevation as the downstream boundary condition. 
As previously noted, the downstream salinity boundary condition for all runs 
was computed from the two profiles presented in Figure 11. The tests for all 
hydrographs except 1980-81, which of course was used to verify the model, 
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should be viewed as estimates of how the current and proposed channels would 
respond to these hydrographs rather than specific simulations of what happened 
during those periods since no attempt was made to represent channel conditions 
at those times, e.g., recall the previous discussion in paragraph 14 concern- 
ing the shoaling of South Pass since the 1973 flood. However, the general in- 
dications should be similar to observations during those periods. Results for 
each hydrograph tested are summarized in subsequent paragraphs. 

47. 19S0-81 hydrograph. The low-flow period experienced began in late 
December 1980 when the river discharge fell below 300,000 cfs on its way to a 
minimum of about 150,000 cfs and continued into mid-February 1981. The re- 
corded Tarbert Landing hydrograph is shown in Figure 31b. The computed time- 
histories of corresponding wedge movement for the 40-, 45-, 50- and 55- ft 
channels are presented in Figure 31a. The maximum intrusion computed during 
the hydrograph reached approximately the position indicated from the steady- 
flow tests from the earlier study for the 150,000-cfs flow for both the 40- 
and 55-ft channel conditions. However, as expected, the time which the wedge 
was computed to be at that location was significantly less than the time which 
the flow fluctuated around 150,000 cfs. As noted In the discussion of model 
calibration, the wedge intrusion for the existing channel is appreciably de- 
layed by the influence of the crossing at mile 47. For the 50- and 55-ft 
channel, the wedge intrusion is primarily influenced by the crossing at 

mile 77. 

48. 1968 hydrograph. Unlike the 1980-81 hydrograph, the 1968 hydro- 
graph experienced a fall and a subsequent rise back to about 300,000 cfs be- 
fore moving to its ultimate low flow of about 160,000 cfs. The recorded dis- 
charges are given in Figure 32b with the computed wedge movement presented in 
Figure 32a. The crossing at mile 47 completely stops the computed wedge for 
the existing and the 45-ft channel. Although this crossing does not apprecia- 
bly retard the movement of the wedge into the system for the 50- and 55-ft 
channel, it has a large Impact on the movement of the wedge out of the system. 
The crossing at mile 63 stops the wedge for the 55-ft channel. 

49* 1955-56 hydrograph. The computed wedge positions corresponding to 
the 1955-56 flows on the Lower Mississippi River are shown in Figure 33a, with 
the discharge hydrograph presented in Figure 33b. With such flows the com- 
puted wedge for both a 50- and a 55-ft channel moves past New Orleans (about 
mile 103) after about 215-220 days. The wedge is stopped by Kenner Crossing 
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Figure 31. Results from 1980-81 hydrograph test 
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Figure 33. Results from 1955-56 hydrograph test 
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at mile 116. As can be seen, the wedge computed with the existing 40-ft chan- 
nel does not match the field data. However, again it should be remembered 
that the channel conditions which currently exist below the Head of Passes are 
quite different from those that existed In 1955-56. Therefore one should not 
expect the computed wedge for existing conditions to match the field data. 
Figure 34 presents wedge profiles along the river for both the 40- and 55-ft 
channels at a time when the concentrations are at their highest point In the 
water column at New Orleans. 
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Figure 34. Wedge Interfaces at the time of minimum 
wedge depth at New Orleans during the 1955-56 

computations 

50. 1953-54 hydrograph. The 1953-54 hydrograph presented in Figure 35b 
resembles the 1980-81 hydrograph but with a lower flow that is maintained for 
a much longer period of time. As shown in Figure 35a for the wedge movement, 
the crossings at miles 47, 63, and 77 have a relatively minor effect on the 
intrusion of the wedge even for the 40-ft channel. Kenner Crossing at 
mile 116 completely stops the wedge for all channels. Figure 36 presents 
wedge profiles along the river at a time when the concentrations are at their 
highest point in the water column at New Orleans. An Interesting observation 
is that the computed wedge for the 55-ft channel comes quite close to matching 
the field data. One might be tempted to conclude that deepening Southwest 
Pass to 55 ft impacts the intrusion of the wedge about as much as If both 
South Pass and Pass a Loutre were restored to their 1953-54 conditions. One 



r? 
IB 



m 

ID 

>- 
o 
o 

CM 



rt 
in 

H 
O 
O 

N 
CM 



in 

> 
G 
2- 



en 

if) 

O 

LLJ 

D 



CO 

to 

O 
Ul 
Q 

*- O 

CM i- 



Z 
< 



2 
< 

—> 

o 

<r> 



« CO 

LU UJ 

-icn 
^cn 



K 

o< 

Q. j£ 

L1J UJ 

DO 

Lurn 

5< 



40 FT CHANNEL (COMPUTED) 



— — *5 FT CHANNEL 

50 FT CHANNEL 

S5 FT CHANNEL 
FIELD DATA 




20 40 60 80 100 120 140 160 180 
TIME, DAYS 

a. Wedge position versus time 



CO 

to 
O 
< 



CO 

m 

O 
D 
< 

co 

CM 





CO 

m 

h- 
O 

o 


co 
w 


co 

in 


co 


m 

IT) 


Q. 
UJ 


1- 
O 


> 
O 


LO 

O 


a 

UJ 


CO 


o 

C\J 


z 


LU 
Q 


o 



r- CM 



t 


T 


in 


m 


7 


Z 


< 


< 


-> 


-J 


o 


o 


^ 


CO 




20 40 60 80 100 1 20 140 1 60 1 80 
TIME, DAYS 
b. Discharge input at Baton Rouge 
Figure 35. Results from 1953-54 hydrograph test 
v ay to substantiate this would be to employ a recently developed version of 
^AEM that has a branching capability. However, since this would constitute 
a substantial modeling effort, arrested wedge computations were made Instead 
in an attempt to provide some insight. 
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Figure 36. Wedge interfaces at the time of 
minimum wedge depth at New Orleans during the 
1953-54 computations 

51. These computations were made by solving a continuity equation for 
the flow in the upper freshwater layer, i.e., 



duh.B 
~djT" 



- Q, 



(26) 



and momentum equations for both the freshwater layer and the stationary bottom 
saline layer, i.e., 



'( h I + h 2) , u du 
dx g dx 



d ' h I + h 2) . u du _ 1 _ 

g dx yhj o 



(27) 



H)S 



p / dx dx yh„ o 



(28) 



with the x-axis drawn in the opposite direction of freshwater flow. 
in the equations are defined as: 

u ■ mean freshwater velocity, fps 
h. - depth of freshwater layer, ft 

B ■ channel width, ft 



Variables 



/, i. 



3 
Q = lateral outflow per unit length, ft /s/ft 

h„ • depth of saline layer, Ft 

2 
g ■ acceleration due to gravity, ft/s 

2 
t = interfacial shear, Ibm/ft/s 

2 2 
Y = mean specific weight for fresh and salt water, lbm/ft /s 

S = bottom slope, ft/ft 

° 3 

Ap = differential density between fresh and salt water, lbm/ft 

3 

p = density of fresh water, lbm/ft 

These equations have been taken from Schijf and Schonfeld (1953). Writing the 
interfacial shear stress as 



Pf i - 
r. - -gi u 2 (29) 



where f , is a nondimensional interfacial friction factor, combining Equa- 
tions 27 and 28 and replacing the derivatives in the equations by centered 
differences results in 



i+1 i h l Bl Q J AX 

hJ +1 B 1+1 hJ +1 B i+1 



i+1 i 
, i+1 . i u + u , i+1 i v 

h = h. „ . , (u - u ) 

I I 2gAp/p 



f.(u +u I Ax 



+ ~ i6 « Ap/p ' Ufx^nf^h* l (31) 



To solve these equations, the freshwater discharge upstream of the wedge must 
be specified and densimetric critical flow is specified at the river mouth, 



v^T 



■ I (32) 
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The solution of Equations 30 and 31 is an iterative one in that u 



1+1 



is ob- 



tained by guessing at h 



i+1 



This value of 



i+1 



is then used to compute a 
value for h, which in turn is used to recompute u . This process con- 



1 



i+1 



tinues until the value of h converges within some specified tolerance, e.g. 
0.01 ft. These computations are patterned after Balloffet and Borah (1985). 
52. Two cases were computed. The first is illustrated below. 
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Here it is assumed that two identical channels with widths of 3,000 ft and 
depths of 40 ft lead to the ocean. Each channel carries 30 percent of the 
total riverflow. As shown below, the second case is for a single 55-ft chan- 
nel leading to the ocean and also carrying 30 percent of the total flow. 
Using values of fAp/o) = 0.02 and f = 0.003 and an assumed water-surface 
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profile representative of the lower Mississippi River, arrested saline wedges 
Were computed for freshwater discharges of 150,000 and 100,000 cfs. At 
70 miles from the ocean end, the freshwater depth was computed to be about 
5 ft greater with the two 40-ft channels when the river discharge was 
150,000 cfs. However, with a discharge of 100,000 cfs the freshwater depth at 
the same location was computed to be only about 1 ft greater with the two 
40-ft channels. The general conclusion then, based upon arrested saline wedge 
theory, is that as the river discharge decreases, the salt wedge created by 
two 40-ft channels approaches the wedge created by a single 55-ft channel. In 
other words, arrested saline wedges created as a result of extremely low riv- 
erflow for channel conditions simulating those before 1973 were approximately 
°f the same strength as those that would be created by a 55-ft channel in 
Southwest Pass today under the same flow conditions. 

53. A comment at this time on how the tip of the wedge is determined as 
the riverflow increases and forces the wedge back toward the Gulf is in order. 
for consistency, when the salinity concentration drops below 9 ppt In a column 
of computational cells the wedge tip Is considered to be located in the adja- 
cent downstream column. Since concentrations on top of the crossings are the 
first to fall below 9 ppt, jumps can be seen in wedge position plots (for 
example, day 148 of the 55-ft channel curve in Figure 35a). This is further 
Illustrated in Figure 37 where plots of the 9-ppt interface on two different 
days are presented. 




90 70 50 30 
MILES ABOVE HEAD OF PASSES 

Figure 37. Determination of the wedge tip 
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54. 1952-53 hydrograph. The minimum flow in the 1952-53 hydrograph 
presented in Figure 38b is about 115,000 cfs and lasted for about 35 days. As 
a result, the computed wedge for all channels (Figure 38a) moved past 
New Orleans to Kenner Crossing, with the wedge from the deepened channels 
remaining there much longer. Plots of the 9-ppt profiles for both the 40- and 
55-ft channels are given in Figure 39 at the time of the minimum depth of the 
wedge interface at New Orleans. 
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Figure 38, Results from 1952-53 hydrograph test 
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Figure 39. Wedge interfaces at the time of mini- 
mum depth at New Orleans during the 1952-53 

computations 

55. 1947-48 hydrograph. The 1947-48 hydrograph contains a relatively 
smooth fall from 300,000 cfs to a low of about 120,000 cfs over a period of 
90 days with a subsequent rise to 220,000 cfs over the next 30 days (Fig- 
ure 40b). As plotted in Figure 40a, the crossing at mile 47 retards the wedge 
in the 40-ft channel for about 25 days with subsequent intrusion to the 
crossing at mile 87 where it remains until the flow increases enough to move 
it back toward the Gulf. However, for the 55-ft channel, the computed wedge 
moves past New Orleans with about a 5- to 10-day retardation period at each of 
the crossings. Wedge profiles along the river are presented in Figure 41. 

56. 1938-39 hydrograph. The 1938-39 hydrograph plotted in Figure 42b 
represents a rapidly varying flow between a high of 300,000 cfs and a low of 
140,000 cfs. The maximum extent of the computed wedge for the 45-, 50- and 
55-ft channels is the crossing at mile 77 (Figure 42a). The major interest in 
this hydrograph is in the rapid manner in which the wedge for both the exist- 
ing and the deepened channels responds to the rapidly varying flow during the 
series of rises and falls that occur. 

57. 1936 hydrograph. The low flow experienced in 1936 extended over a 
rather lengthy period with flow dropping below 300,000 cfs in early June on 
its way to a minimum flow of 91,600 cfs some 95 days later (Figure 43b). The 
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b. Discharge input at Baton Rouge 
Figure 40. Results from 1947-48 hydrograph test 
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Figure 41. Wedge interfaces at the time of minimum 
wedge depth at New Orleans during the 1947-48 

computations 
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b. Discharge input at Baton Rouge 
Figure 42. Results from 1938-39 hydrograph test 
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b. Discharge input at Baton Rouge 
Figure 43. Results from 1936 hydrograph test 

computed wedge movements in response to the 1936 hydrograph are plotted in 

Figure 43a. The maximum intrusion for all channels is Kenner Crossing. The 
effect of the crossing at mile 47 on the wedge is clearly evident. As pre- 
viously discussed in paragraph 53, Figure 37 illustrates the manner in which 
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the tip of the wedge is determined. The primary wedge for the 1936 hydrograph 
moves out rapidly but pockets of high-salinity water are temporarily trapped 
on the upstream side of channel crossings. Figure 44 presents a plot of the 
9-ppt profiles for both the 40 and 55-ft channels. 
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Figure 44, Wedge interfaces at the time of 
minimum wedge depth at New Orleans during 
the 1936 computations 

Duration Results 



58. These hydrograph tests provide information on the extent of the 
wedge intrusion for the 40- , 45-, 50-, and 55-ft channels under a wide range 
of flow hydrographs that have occurred in the past. In addition, results from 
the tests can be used to estimate the duration of wedge intrusion, based upon 
a concentration of 9 ppt, beyond Port Sulphur (mile 40), Algiers (mile 96), 
and New Orleans (mile 103) , Increases in duration as a result of the deepened 
channels are presented in Tables 2-4. Figures 45 and 46 contain plots for 
both the 40- and 55-ft channels of duration of wedge intrusion versus river 
miles for the 1953-54 and 1936 hydrographs, respectively. The influence of 
the various river crossings is responsible for the sharp changes in slope on 
some of the curves. For the 55-ft channel, the 1953-54 hydrograph results in 
a wedge that stays past New Orleans for 67 days (Figure 35a or Figure 45) , 
The conditions at New Orleans generated by the 1952-53 (Figure 38a) and the 
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1936 (Figure 43a or Figure 46) hydrographs for a 55-ft channel appear to be of 
approximately equal severity, with the wedge remaining past New Orleans for 
about 40 days. Hydrographs such as the 1980-81, 1968, and 1938-39 would not 
create problems for locations beyond the crossing at mile 77. 
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PART VII: SILL TESTS 

59. Results from the hydrograph tests reveal that deepened channels up 
to 55 ft will significantly increase the extent and duration of salinity in- 
trusion. As a possible means of either stopping or at least severely retard- 
ing the intrusion of the salt wedge during critically low-flow periods, a sill 
constructed from bed material that is primarily fine sand has been proposed at 
the crossing located at about mile 63. The sill would have a triangular shape 
with IV on 40H side slopes and a maximum top elevation of -55 ft. The stabil- 
ity of the increased height of the crossing was determined using the HEC-6 
computer program "Scour and Deposition in Rivers and Reservoirs" (HEC 1977); 
whereas its impact on the extent of salinity intrusion was of course deter- 
mined by employing the flow-salinity model, LAEM. The sill stability tests 
were conducted by Copeland (1983), and the following discussion was taken from 
his summary of the results of the computer runs. 

Sill Stability 

60. The study reach extended from Venice, at river mile 10.7, to 

New Orleans, at river mile 102.7. An HEC-6 model set up by LMN was used as a 
base for geometric and bed material data. Additional cross sections were 
taken from the 1973-75 hydrographic survey. The Toffaleti sediment transport 
function was used in the model because it favorably reproduced sediment mea- 
surements from Tarbert Landing (river mile 316). Sediment input into the nu- 
merical model was initially determined from calculated hydraulic parameters. 
The sediment input was adjusted to maintain an essentially stable bed eleva- 
tion at the upstream cross section. Cross sections identified as erosion- 
resistant were not allowed to erode in the numerical model. A report by Kolb 
(1963) was used to make these determinations. The downstream water-surface 
fating curve was determined from daily stages at Venice and discharges routed 
from Tarbert Landing and Carrol lton. Hydrographs were developed from daily 
discharge records at Tarbert Landing and were limited to a maximum of 
1,250,000 cfs; this accounted for operation of the Bonne Carre and Morganza 
^loodways. The Pointe a la Hache Relief Outlet, located between river miles 
33 and 44, allows some flow to leave the river when the water-surface eleva- 
tion exceeds 7 ft. At a maximum discharge of 1,250,000 cfs, LMN estimated 

» 
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that the outflow would only be about 30,000 cfs or 2 percent of the total. 
The effect of the outlet was considered insignificant and was therefore not 
simulated in the study. 

61. The numerical flow model was verified by matching recorded stages 
at eight gages for a range of steady-state discharges. This was accomplished 
by adjusting Manning's roughness coefficients. Adopted values varied with 
discharge and ranged between 0.013 and 0.028 from Venice to river mile 50 and 
between 0.020 and 0.035 from river mile 50 to New Orleans. Calculated water- 
surface elevations compared favorably with measured stages as shown in 
Figure 47. 

INNER HARBOR 
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GAGE READINGS 
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Figure 47. Water-surface profiles from HEC-6 

model 

62. The sediment model was verified by comparing measured and calcu- 
lated sediment transport rates and bed changes. Calculated transport rates at 
river mile 76 compared favorably with measured values at the Belle Chase gage 
(Figure 48), Data from the 1963-65, 1973-75, and 1982 hydrographic surveys 
indicated that bed elevations in the study reach did not change significantly. 
Eight years of record (1975-1982) were simulated in the model; bed changes in 
the vicinity of the sill were insignificant (Figure 49). The model was not 
sufficiently verified below river mile 50 where significant scour and 
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deposition were calculated, but these results did not affect the calculations 
in the area of interest. 

63. The verified numerical model was used to determine the stability of 
the sill. Steady-state discharges ranging between 200,000 and 800,000 cfs 
were simulated for a period of 150 days; the following depths were eroded off 
the top of the sill: 



Discharge, cfs 

200,000 

300,000 
400,000 
450,000 



Erosion, ft 

0.02 

1.4 

2.0 

5.0 



Rates of erosion for higher discharges are shown in Figure 50. At a discharge 
of 800,000 cfs, the sill washed out in about 5 days. The sill may be consid- 
ered stable when the discharge is less than 400,000 cfs. 
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Figure 50. Erosion of triangular sill (IV on 40H) 

at mile 63.4 

Impact of Sill on Salinity Intrusion 



64. Of all the hydrographs tested, the 1953-54 hydrograph produced the 
most severe salinity conditions at New Orleans, although the low flow of rec- 
ord occurred in 1936 and historical records imply that salinity conditions 
were more severe in 1936. However, again it should be remembered that channel 
conditions in 1936 and 1953 were probably not the same, whereas, the same 
channel was used in the model for all runs. In addition, values for the re- 
corded river flows employed could be in error by as much as perhaps 10-15 per- 
cent. Therefore the 1953-54 flow conditions at Baton Rouge, as well as the 
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1936 conditions, were used to test the effect of various sills at mile 63 on 
salinity intrusion. The wedge intrusion for a sill height of -55 ft is pre- 
sented in Figures 51 and 52 for the 1953-54 conditions and 50- and 55-ft chan- 
nels, respectively. A similar plot for a 45-ft channel but with the sill at 
-60 ft is presented in Figure 53. Likewise, the wedge intrusion for a sill 
height of -55 ft is presented in Figures 54 and 55 for the 1936 conditions and 
50- and 55-ft channels, respectively. A similar plot for a 45-ft channel and 
sills at heights of -65, -60, and -55 ft is presented in Figure 56. Although 
a sill does not completely stop the wedge, the movement of the wedge is re- 
tarded such that the overall intrusion with the -55-ft sill in a 55-ft channel 
closely follows that computed with the 40-ft channel. Decreases in duration 
beyond Port Sulfur, Algiers, and New Orleans, as a result of a -55-ft sill in 
a 55-ft channel, are presented in Table 5. 

Effect of Sill Discretization on Salinity Intrusion 

65. As discussed in paragraph IS, the earliest runs of LAEM utilized a 
spatial step of 1 mile, whereas the runs discussed here employed a Ax of 
2 miles in the geometric representation of the system. Therefore a sill con- 
structed on the natural crossing at mile 63 is assumed to be 2 miles wide. 
The question of the effect of the width of the sill on salinity intrusion has 
been posed. Based upon the fact that essentially no differences in the com- 
putations were observed between runs in which river crossings were represented 
by 1 mile and then 2 miles, it is not believed that the width of the sill has 
an appreciable effect on the flow-salinity fields. 
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Figure 51. Effect of sill at mile 63.4 on salt wedge 
intrusion in the 50-ft channel for 1953-54 hydrograph 
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Figure 52. Effect of sill at mile 63.4 on salt wedge 
intrusion in the 55-ft channel for 1953-54 hydrograph 
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intrusion in the 45-ft channel for 1953-54 hydrograph 
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Figure 55. Effect of sill at mile 63.4 on salt wedge 
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PART VIII: SUMMARY AND CONCLUSIONS 
Summary 

66. Since there is a concern that deepening Southwest Pass could sig- 
nificantly influence the salinity intrusion in the Lower Mississippi River, 
which in turn could have a serious impact on the supply of fresh water to 
New Orleans and other localities along the river, a dynamic analysis of the 
salinity intrusion phenomena using actual channel cross sections was autho- 
rized by LMN. To accomplish this, a numerical model called LAEM, which is a 
laterally averaged model that couples the computation of the flow and salinity 
fields through the influence of the salinity on the water density, has been 
employed. 

67. After verifying the model using 1980-81 field data, several his- 
torical hydrographs furnished by LMN were run and results for 40-, 45-, 50-, 
and 55-ft channels have been presented. In addition, results from a series of 
steady-flow runs, obtained by using the model as previously verified without 
several of the river crossings below New Orleans, are presented and tend to 
support the intrusion curves developed several years ago by LMN. 

68. The final application of the model was to test the impact on the 
intrusion of the salt wedge of increasing the height of the natural river 
crossing at about mile 63 by dredging the upstream channel and depositing the 
material on the crossing. An assessment of the stability of a sill built from 
natural river sediment was accomplished through an application of the HEC-6 
computer program. 

Conclusions 

69. In general, it is concluded that the numerical model, LAEM, pro- 
vides a very useful tool for assessing the Impact of major changes in channel 
geometry on the intrusion of salinity. Verification of LAEM for the Lower 
Mississippi River produced a model whose computed results were in good agree- 
ment with the limited field data available and responded well to the various 
boundary data imposed. Particular conclusions from the salinity-intrusion 
computations reported herein follow: 

a. Only the 1980-81 field data could be utilized to verify the 
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model since the shoaling of South Pass since 1973 has changed 
the manner in which the salt wedge is formed by saline water 
from the Gulf. 

b. Computed wedge intrusion for earlier steady-flow tests In both 
the 40— and 55-ft channels was generally supportive of the 
intrusion curves developed by LMN. 

c. Hydrograph tests for the 1930-81, 1968, 1955-56, 1953-54, 
1952-53, 1947-48, 1938-39, and 1936 low-flow periods provide 
useful information on duration of wedge intrusion for these 
hydrographs and insight concerning the importance of hydrograph 
shape and channel crossings on the extent and duration of wedge 
intrusion. 

d. Channel crossings, particularly the ones at miles 47, 77, 87, 
and 116, significantly limit the extent and duration of salin- 
ity intrusion in the Lower Mississippi River. 

e. Results from runs with deepened channels revealed that deepen- 
ing the channel will significantly increase the extent and 
duration of salinity intrusion. 

f. The 1953-54 hydrograph resulted in the greatest intrusion and 
duration of the computed salt wedge, although the low flow of 
record on the Lower Mississippi River occurred during 1936 and 
historical data imply that more severe salinity conditions ex- 
isted in 1936. 

g. An interesting observation from the computations for the 
1953-54 flows is that the wedge position for the 55-ft channel 
followed quite closely the field data collected during that 
time. Therefore, with channel conditions below the Head of 
Passes as they currently exist, the impact on salinity Intru- 
sion of deepening Southwest Pass to 55 ft may be to essentially 
return the wedge movement back to pre-1973 conditions before 
South Pass was allowed to shoal. Computations from an arrested 
wedge model tended to substantiate this conclusion. 

h. Man-induced alterations to create greater heights at river 

crossings, i.e., to create a sill in the river, appear to be an 
effective means of limiting wedge intrusion during critically 
low-flow periods. For 1953-54 flow conditions and a 55-ft 
channel, a sill at mile 63 with a top elevation of -55 ft would 
result in salinity conditions at New Orleans similar to those 
that occur with the existing 40-ft channel. 

70. Results from the sill stability tests with a -55-ft sill in a 55-ft 
channel indicate that the deposited material on the crossing will be eroded 
fairly quickly as the riverflow Increases above 400,000 cfs, e.g., a constant 
riverflow of 700,000 cfs for 10 days will erode about 9 ft of the sill height. 
For rlverflows below 400,000 cfs, the sill will be stable. These results 
imply that the sill will have to be reconstructed before each critical low- 
flow event. 
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Variables Used 



Table I 

in Computation of t. 

















River Mile 
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dh 2 /dx 
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Table 2 
Increase of Duration of Wedge Intrusion (9 ppt), with a 45-ft C h ann el, 

Beyond River Locations, Days 
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Table 3 
Increase of Duration of Wedge Intrusion (9 ppt), with a 50-ft Channel, 

Beyond River Locations, Days 



























Hydrograph 


Port 


Sulfur, 


Mile 


40 


Algiers, Mile 


96 


New 


Orleans , 


Mile 103 


1980-81 




7 






















1968 




29 






















1955-56 




30 








— 








— 




1953-54 




8 








19 








20 




1952-53 




28 








22 








22 




1947-48 




29 






















1938-39 




28 






















1936 




26 








9 








9 





Table 4 
Increase of Duration of Wedge Intrusion (9 ppt), with a 55-ft Channel, 

Beyond River Locations, Days 

























Hydrograph 


Port 


Sulfur, 


Mile 


40 


Alg 


ie 


rs, Mile 


96 


New Orl 


eans, Mile 103 


1980-81 




14 




















1968 




46 




















1955-56 




53 










— 






— 


1953-54 




10 










40 






27 


1952-53 




49 










30 






31 


1947-48 




35 










14 






12 


1938-39 




43 




















1936 




30 










14 






14 



Table 5 

Decrease of Wedge Intrusion (9 ppt), in a 55-ft Channel, with a Sill at 

Elevation -55 at Mile 63, Beyond Particular Locations, Days 

Hydrograph" Port Sulfur, Mile 40 Algiers, Mile 96 New Orleans, Mile 103 
1953-54 30 31 

1936 17 18 



